Assessment of soil heavy metal pollution and associated ecological risk of agriculture dominated mid-channel bars in a subtropical river basin

The elevated concentrations of heavy metals in soil considerably threaten ecological and human health. To this end, the present study assesses metals pollution and its threat to ecology from the mid-channel bar’s (char) agricultural soil in the Damodar River basin, India. For this, the contamination factor (CF), enrichment factor (EF), geoaccumulation index (Igeo), pollution index, and ecological risk index (RI) were measured on 60 soil samples at 30 stations (2 from each station, i.e., surface and sub-surface) in different parts of the mid-channel bar. The CF and EF indicate that both levels of char soil have low contamination and hence portray a higher potential for future enrichment by heavy metals. Moreover, Igeo portrays that soil samples are uncontaminated to moderately contaminated. Further, pollution indices indicate that all the samples (both levels) are unpolluted with a mean of 0.062 for surface soils and 0.048 for sub-surface soils. Both levels of the char have a low potentiality for ecological risk with an average RI of 0.20 for the surface soils and 0.19 for the sub-surface soils. Moreover, Technique for order preference by similarity to ideal solution (TOPSIS) indicates that the sub-surface soils have lower pollution than the surface soils. The geostatistical modeling reveals that the simple kriging technique was estimated as the most appropriate interpolation model. The present investigation exhibits that reduced heavy metal pollution is due to the sandy nature of soils and frequent flooding. However, the limited pollution is revealed due to the intensive agricultural practices on riverine chars. Therefore, this would be helpful to regional planners, agricultural engineers, and stakeholders in a basin area.

The tropical rain-fed river Damodar and its basin area experience the southwest monsoonal regime 43 . The mean annual temperature in the basin area lies between 18° and 24 °C in winter and 29°-35 °C in summer with a mean annual rainfall of 100-200 cm. Based on the census of India 49 , about 17.25 million people live with a 738 km 2 population density in the basin area. Major cities like Dhanbad, Bokaro, Durgapur, Asansol, and Bardhaman are located in the DRB 41 . Intensive agriculture is practiced in the basin area including its chars. The gross irrigation command area of this basin is 5690 km 2 out of which 300 km 2 of land in the upper basin is being irrigated and 3640 km 2 of irrigation potential land is created by DVC dams in the basin area 50 . National Waterway 35, from Krishak setu (Bardhaman on State Highway No. 8) to the confluence with Hooghly River near Purba Basudebpur, 46% of land along the Damodar River is characterized as agricultural, and on average 58% of Bardhaman district total population belongs to the agricultural population. The agricultural effluents with numerous pollutants are discharged into the river at different points from the agricultural field of its course directly or indirectly, leading to an ecological risk for the river system. Sampling design and data collection. This study has been conducted based on the physico-chemical parameters of soils that have been collected from the mid-channel bars agricultural field of Damodar River. A total of 60 soils have been collected from the selected mid-channel bars (chars) agricultural at 30 different locations in 2021 (Fig. 1). Two soil samples were collected from each station i.e. one is from the surface level and another is from the sub-surface. As the char soils (fluvisols) are annually flood inundated, they do not exhibit mature soil profiles. Hence, soil profile depth has been considered as the criteria to segregate the surface and sub-surface soils. The depth of collected surface soils is 0-0.1 m and the sub-surface is 0.7-1 m (Table S1). Before sampling, the uppermost layer of soil has been removed for top-level soil collection, and for sub-surface soil collection, a digging bar (Sabal) is used to dig and collect soil. A quadrant pattern of the grid is used to collect soil samples in the study area 51 . Based on this method, a 1 m 2 grid has been considered for taking soil samples and four samples from the surface soil in each location were collected and mixed well to prepare one soil surface sample. A similar approach was adopted for preparing the sub-surface soil samples. From each site, a well-mixed 500-600 g soil has been collected in an air-tight plastic polythene bag or pouch. After that, all the collected soils were corked, labeled carefully, and packed in another airtight bag. After collecting all the soils from the field, they were given for testing in the laboratory of Vivekananda Institute of Bio-Technology (Soil, Water and Manure Testing Laboratory) at Sri Ramkrishna Ashram in Nimpith, South 24-Parganas, West Bengal. The physico-chemical parameters such as pH, EC, OC, CEC, soil texture (sand, silt, and clay), soil moisture, and HM like Cu, Zn, Mn, Fe, and S were tested for each soil sample. We collected soils from char in pre-monsoon (February) because the concentration of pollution can be detected easily and efficiently. In other seasons (monsoon and post-monsoon), the concentration of pollution is low due to the monsoon and higher precipitation 52,53 . Moreover, we restricted our study only to pre-monsoon season because some parts of the char (mostly the northwestern part) are low-lying and often become submerged under monsoon flood water.
Pedo-chemical assessment and estimation of heavy metal pollution. Pedo-chemical assessment has been executed using a robust methodology starting from the thematization, development of the research design for the collection of the data, and representation of the data with suitable geospatial and geo-statistical methods (Fig. 2).
Heavy metal pollution indices and ecological risk estimation. Several heavy metals contamination and ecological risk assessing methods such as contamination factor (CF), contamination degree (Cd), potential contamination index (Cp), enrichment factor (EF), geoaccumulation index (I geo ), pollution index (PI), pollution load index (PLI), ecological risk index (RI), and potential ecological risk coefficient (Er) have been used for measuring HM pollution in soil and their risk to ecology from the studied soils samples. The CF, Cd, and Cp are used to assess HM contamination in soil 33,34,54 while EF is used to measure HM enrichment in soil 55 and to find out the potential sources of HM as the human or geogenic origin 10,33 . Similarly, I geo is used to judge the level of metal accumulation in soil 10 . Moreover, PI and PLI are used to evaluate the pollution level in soil in a specific zone or site 33,56,57 . Finally, the ecological risk from HM is assessed through Er, RI and MRI to judge the sensitivity of different communities 30,54 . These methods employed in the present investigation are widely used for evaluating soil pollution and ecological risk from the PC parameters of soil. The details of these methods and algorithms are given in Table 1.
TOPSIS MCDM model. Multi-criteria decision-making (MCDM) is a significant method to determine a problem that has multiple conflicting criteria in decision-making 58,59 . TOPSIS (Technique for Order of Preference by Similarity to Ideal Solution), an MCDM technique was invented by Hwang and Yoon in 1981. To calculate the optimum alternative, it accounts for both the shortest distance from the ideal solution and the longest distance from the anti-ideal solution at the same time 60 . Therefore, it is a useful technique to solve decision-making problems in the real world and is a widely used MCDM technique for prioritization purpose 61 . Hence, in the present investigation, the TOPSIS has been implemented to evaluate the priority of surface and sub-surface soil pollution. To avoid subjectivity, the weights of the variables used in the TOPSIS method, have been defined based on The soil pollution and ecological risk assessing indices such as Cd, Cp, EF, I geo , MPI, PLI, Er, RI and MRI have been analyzed by TOPSIS for prioritizing the locations of the surface and sub-surface soil pollution for better management practices. The surface and sub-surface soil of the mid-channel bar is selected as the alternatives, and 16 variables related to soil pollution and ecological risk assessment are used as the criteria in the matrix table (Fig. 3).
Geostatistical analysis. The geostatistical model is used to analyze and predict the values related to spatial or spatiotemporal phenomena 63 . Geostatistical modeling has been used for the present study for error minimizing and finding the best-fit model to show the spatial distribution of variables. The best-fit model depends on the data nature and its spatial or spatiotemporal distribution. We applied geostatistical modeling for the spatial distribution of PC parameters and index values of soils. The interpolation methods like simple kriging (SK), ordinary kriging (OK), and inverse distance weighting (IDW) have been used for assessing the spatial distribution of PC parameters and index values of soils of the studied chars. The working principle and algorithms of SK, OK and IDW and best-fit model assessing tools like semi-variogram 18 coupled with other statistical measures such as mean square error (MSE), mean error (ME), root mean square error (RMSE), root mean square standardized error (RMSSE), and average standard error (ASR) used in the present study are provided in supplementary (ST1). We used ArcGIS software (v. 10.2) for running the interpolation models and testing the best-fit ones.

Results
Pedo-chemical characterization. Soil textural distributions. Textural analysis of soil is an important factor for studying the PC or HM pollution in soil because the concentration of PC elements, essential nutrients, and other living organisms depends on soil texture. The textural analysis exhibits that the majority of soils are mostly sandy and classified as fluvisols, a type of entisols developed by the fluvial deposits following annual flood events. In general, the sizes of soil particles affect the concentration of HM. The sizes of soil particles decrease with an increase in the concentration of HM in soil 64 Table 2). Thus sub-surface soils are sandier than the surface. Furthermore, Shepard Triangle Diagram 65 is used to classify the soils (Fig. 4a,b). This diagram showed that 83.34, 13.33 and 3.33% of soils were sandy loam, sandy clay loam and loam of surface level. Besides, 90, 3.33 and 6.67% of soils were sandy loam, sandy clay loam and loam of sub-surface level respectively.
Soil pH and EC. The pH is an important chemical property of soil as it controls precipitation and adsorption which are the principal mechanisms of metal retention and affect metal dynamics in soil 64 . It plays a vital role in soil fertility as it controls the movement of HM and activates the macronutrients and trace elements in soil 26,66 . The present study found that surface soils were moderately acidic while the sub-surface was slightly acidic (surface pH = 5.49 and sub-surface pH = 6.11). The mean pH values of surface and sub-surface soils indicate that surface soil has nutrient deficiencies and toxicities (iron, manganese, and aluminum) and other elements (calcium, magnesium, nitrogen, phosphorous, and potassium) become less available for absorption by plants whereas subsurface soil is an ideal soil condition for most plants. The pH of surface soils ranges from 4.48 to 6.43 while in the sub-surface soils, it ranges from 4.54 to 7.01 (Table 2). Based Ficklin-Caboi diagram, 76.67% of surface soil and 16.67% of sub-surface soils lie under acid high metal while 23.33% and 83.33% lie under near natural high metal respectively (Fig. 5a, b). Electrical conductivity (EC) is used to indicate the salinity of soil and water 67,68 . The present study found a low concentration of EC in the soils (0.05-0.69 dsm −1 with a mean of 0.20 dsm −1 for the surface soils and 0.03-0.37 dsm −1 with a mean value of 0.11 dsm −1 for the sub-surface soils ( Table 2 and Fig. 6a,b). This result reveals that the surface soils are more saline than the sub-surface soils.
Organic carbon. Organic carbon (OC) is also a significant parameter of agricultural soil because it maintains the bio-availability of HM in soil 29 . A low concentration of OC in soil reduces the soil microbial diversity and www.nature.com/scientificreports/ Table 1. Heavy metal pollution indices and ecological risk estimation.

Method Equation Description
Scale limit and classification system References C o -observed metal concentration in soil; C b -background value CF < 1 denotes low, 1 ≤ CF < 3 for moderate, 3 ≤ CF < 6 for considerable and 6 ≤ CF implies very high 10,26,27 Degree of contamination Cd = n i=1 CF CF-Contamination factor for metal; n-Number of HM  30,33 Enrichment factor EF = (Cx/Fe)sample (Cx/Fe)background C x sample and C x background -concentration of metal in soil and background environment; Fe sample and Fe background -concentration of Fe in soil and the background environment EF ≤ 1 for no enrichment, 1 < EF ≤ 3 for minor , 3 < EF ≤ 5 for moderate, 5 < EF ≤ 10 for moderately severe, 10 < EF ≤ 25 for severe enrichment, 25 < EF ≤ 50 for very severe and EF > 50 for extremely severe [28][29][30] Geoaccumulation index (I geo ) I geo = log 2 Cn 1.5×Bn C n -metal concentration in soil sample; B n -background value I geo ≤ 0 represents practically uncontaminated, 0 ≤ I geo < 1 for uncontaminated to moderately contaminated, 1 ≤ I geo < 2 for moderately contaminated, 2 ≤ I geo < 3 for moderately to heavily contaminated, 3 ≤ I geo < 4 for heavily contaminated, 4 ≤ I geo < 5 for heavily to extremely contaminated and 5 ≤ I geo for extremely contaminated [26][27][28] Pollution index PI = CF average -average Contamination factor of metals; CF maximum-maximum contamination factor of metals PI < 0.7 implies unpolluted, 0.7 < PI < 1 for slightly polluted, 1 < PI < 2 for modest polluted, 2 < PI < 3 for heavily polluted and PI > 3 for severely polluted 32 Modified pollution index MPI = EF average -average Enrichment factor of metals; CF maximum -maximum Enrichment factor of metals MPI < 1 indicates unpolluted, 1 < MPI < 2 for slightly polluted, 2 < MPI < 3 for modest polluted, 3 < MPI < 5 for modest heavily polluted, 5 < MPI < 10 for heavily polluted and 10 < MPI for severely polluted 32 Pollution load index  Table 2 and Fig. 6a,b). These results indicate that the microbial diversity and biomass of sub-surface soils are lower than those of surface-level.
Soil moisture. Soil moisture also impacts the soil organisms and interacts with contaminants of soil. Hence, it affects the solubility of HM and the bioavailability of soil. It was found that the availability of soil moisture in surface soils ranged from 0.80 to 11.40% with a mean of 2.88% whereas it varied from 0.60 to 12% with a mean of 2.31% in the sub-surface soils (Table 2). Thus, the percentage of soil moisture concentration in the surface soils is more than the sub-surface level (Fig. 6a,b).
Cation exchange capacity. Soil CEC is a significant chemical property that reflects soil functions like structural stability, nutrient accessibility, pH, and reaction to fertilizers 64 . It is a measure of the soil's capability to bind interchangeable cations 64,66 . Our study findings indicate that the CEC (meq/100gm) of surface soils was found to be 2-5.44 with a mean of 3.99 while it was found to be 2.52-11.48 with a mean of 6.77 in the sub-surface soils ( Table 2 and Fig. 6a,b). These results also indicate that the CEC of sub-surface soils has more ability to bind or hold the exchangeable cations than surface soils.  (Fig. 6a,b).  (Fig. 7a,b).
Contamination degree (Cd), and potential contamination index (Cp). The spatial distribution of Cd of surface soils was found to be 0.095-0.188 with a mean value of 0.134 whereas it was found to be 0.081-0.174 with a mean value of 0.114 in the sub-surface soils (Fig. 7a,b). The Cd of the studied HM is very low in all the soils. The study exhibits that all the Cd values from both level soils lie below 8C degrees which indicates that the metals have a   (Fig. 7a,b). Based on the calculated EF mean values, it is observed that the EF values of subsurface soils' HM are higher than the surface soils. It is also observed that the EF value of Mn (in both levels) is high among the other investigated metals. According to their EF mean values, the investigated Cu, Zn, and Mn HM of surface soils are classified as severe, very severe, and extremely severe respectively. However, at the subsurface level, they are classified as very severe, very severe, and extremely severe (Table 2). Moreover, they are ranked in the order of Mn > Zn > Cu for both levels of soil.
Geoaccumulation index (I geo ). The study found that the I geo value of surface soils Cu, Zn, Fe, and Mn was 0.0020-0.0069, 0.0027-0.00114, 0.00011-0.00033 and 0.010-0.025 respectively whereas they were 0.0025-0.0065, 0.0032-0.00162, 0.00009-0.00022 and 0.007-0.020 in the sub-surface soils in the study. Therefore, these results show that all the calculated I geo values lie below 1 (< 1) and Class 1 which indicates that soils are uncontaminated to moderately contaminated. The calculated I geo mean values of HM reveal that the I geo values of Zn, Fe, and Mn of sub-surface soils are lower than the sub-surface level, although the I geo value of Cu is higher than  www.nature.com/scientificreports/ the surface. Moreover, based on the mean I geo values of HM, they were ranked in the order of Mn > Zn > Cu > Fe in both levels of soils (Fig. 8a,b).

Pollution index and pollution load index.
In this study, the calculated pollution index values of surface soils ranged from 0. 040 to 0.093 with a mean value of 0.062 while it ranged from 0.031 to 0.073 in the sub-surface soils (Table S2 and Fig. 8a,b). Thus, based on these results, the pollution index values of surface soils are higher than the sub-surface. It indicates that the pollution level is also high in surface soils. Moreover, according to the classification of pollution index values, all the soils' pollution index indices lie below 0.70 indicating the unpolluted nature of all the samples (from both levels). The PLI was found to be 0.11-0.020 in surface soils while 0.009-0.17 in the sub-surface soils (Fig. 8a,b). Hence, these results indicate that the pollution load of HM in surface soils is more than in the sub-surface.
Ecological risk assessment. The ecological risk index (RI) shows that the RI values of surface soils ranged from 0.15 to 0.26 while it ranged from 0.14 to 0.25 in the sub-surface soils (Table S3). The RI values of surface soils are higher than the sub-surface. Hence, the ecological risk from the surface soil is higher than the subsurface. Moreover, all the soils' RI values lie under 150 indicating that all the samples from both levels have low ecological risk (Fig. 8a,b). The E i r values of Cu, Zn, and Mn in the surface soils were 0.050-0.172, 0.014-0.057, and 0.051-0.124 respectively. However, they were found to be 0.061-0.161, 0.016-0.081 and 0.036-0.098 in the sub-surface soils (Table S3). The mean values of E i r show that the E i r values of Zn and Mn of surface soils are higher than the sub-surface but Cu is lower than the sub-surface soils. It is also found that E i r value of Cu (in both levels) is high among the other investigated metals. The mean E i r values of Cu, Zn, and Mn, E i r values of each HM in all  www.nature.com/scientificreports/ the soils from both levels were less than 40 and classified as low potential ecological risk. Therefore, the studied HM in soils has a low potentiality for ecological risk. Moreover, they are ranked as Cu > Zn > Mn for both levels of soil (Fig. 9).

Prioritization of surface and sub-surface soil using TOPSIS. Prioritization of surface and sub-sur-
face soil is vital for the management and mitigation of soil pollution and its ecological risk using TOPSIS 14 . Pollution and ecological risk assessing indices such as Cd, Cp, EF, I geo , MPI, PLI, Er, RI and MRI have been analyzed using the TOPSIS MCDM technique for evaluating the priority of surface and sub-surface soil pollution. The performance rank (Pi) portrays that a higher value of Pi indicates lower soil pollution and vice versa. Hence, rank 1 is the best location while rank 30 is the worst location in terms of soil pollution. The calculated Pi values of surface and sub-surface soils are mentioned in supplementary (Table S4). The final relative closeness score of the ideal solution was found to be 0.29-0.81 with a 0.65 mean value in the surface soils while it was found to be 0.04-0.99 with a 0.94 mean value in the sub-surface soils. These results indicate that the sub-surface soils have lower soil pollution than the surface soils. Furthermore, the final score or Pi values are classified into five groups such as very low (< 0.4), low (0.4-0.5), medium (0.5-0.6), high (0.6-0.7), and very high (> 0.7) (Fig. 10a,b).
Geostatistical modeling. This study found that all the semi-variogram models (spherical, exponential, Gaussian, and circular) have arisen as best fit models for OK and SK interpolation methods for some elements as there was no difference in the ME, MSE, RMSE, ASR, and RMSSE values (  www.nature.com/scientificreports/ observed as the best-fit model for both levels of soils but for sand and clay, SK was observed for surface soils. And OK for sub-surface soils. However, in case of Fe, inverse results were found. For moisture, the spherical semivariogram model of the SK interpolation technique is the best fit for both levels of soil. Further, for EC and Silt, the circular and exponential/Gaussian semi-variogram models of the OK interpolation technique are the best-fit model for surface soils while Gaussian of SK and OK is the best-fit model for surface soils respectively. Hence, based on these results, the kriging interpolation technique (OK and SK) with the studied semi-variogram model provides better performance for each variable (Table 3). Furthermore, 61.54% of surface soil sample variables and 76.92% of sub-surface soil sample variables fit with the SK interpolation technique and the rest with the OK interpolation technique. So, the SK technique is expected as the most accurate interpolation model for OC, CEC, sand, clay, moisture, Cu, Mn, and B in the surface soils while it is for pH, EC OC, CEC, moisture, Cu, Zn, Fe, Mn and B for sub-surface soils as compared to the other studied techniques (OK and IDW). The elements like a nugget, range, partial sill, and lag size value of best-fit semivariogram models are extracted by using ArcGIS (version 10.2). Later, sill, nugget/sill, and the effect of nugget/sill are also calculated from obtained values for this study (Table 3).

Discussion
The present study assessed the PC parameters of collected surface and sub-surface soils of mid-channel bar's agricultural fields in the Damodar River for measuring HM pollution and their ecological risks. This study observed that the HM concentration was low as per the permissible limit of HM based on Salem et al. 66 . Moreover, it was also observed that the Zn, Fe, and Mn concentrations in surface soils are relatively low than in the sub-surface soil. The fractionation of the HM may have been influenced by pH, OC, CEC, and moisture which all depend on soil structure 69,70 . This study observed that the surface soils ranged from strongly acidic to slightly acidic and sub-surface soils varied from strongly acidic to neutral. Although, Rani et al. 29 found the alkaline nature of agricultural soil in Delhi, India and they mentioned that low pH has high solubility of HM in soil. Lower concentrations of pH in soil may decrease HM adsorption and increase soil mobility 70 . Moreover, the study suggests most of the surface soils are contaminated with acidic high metals while most of the sub-surface soils are near natural high metals. Kumar and Singh 10 found neutrals or alkaline with high metals types of sediment in a tropical Ajay River basin system in India, which is closely associated with the present study findings. The major causes of soil acidity are the excess leaching of nitrate, using nitrogen-based fertilizers in the soil, and removing plant and animal products 71,72 . The leaching process accelerates as a result of the sandy soil. Acidic soil impacts agricultural productivity and sustainable farming system by decreasing the availability of essential nutrients and increasing the toxicity of trace elements 73 . The low concentration of OC in the soils was also observed. In addition, OC and the percentage of soil moisture concentration in surface soils were higher but the CEC was lower than in the sub-surface soils. The CEC in soils indicates that sub-surface soil can hold the exchangeable cations more than surface soil. Yu et al. 70 mentioned that reducing soil CEC may increase bioavailability in soil. Moreover, due to the sandy soil texture, the concentrations of these physico-chemical parameters in the midchannel bar soils may be found low as these PC or HM are washed out with run-off as agricultural effluents into the river Damodar (Fig. 11a). These mixing agricultural effluents increase the pollution level of Damodar River water due to the growth of intensive agriculture on the chars (Fig. 11b,c) and the threat to the riverine ecology as well as the riverine biota to some extent. Our study findings are supported by Sarkar et al. 52 .
The analysis of HM contamination shows a low contamination level for CF, Cd, and Cp in both the surface and sub-surface soil. Perumal et al. 33 reported a similar finding in their study on the Thondi coast, Palk Bay, South India. Further, the geoaccumulation index indicates that all the soils from both levels are uncontaminated to moderately contaminated by HM. Our study findings are also supported by Sarma Bora et al. 74 . Although Perumal et al. 33 found that Fe and Mn were extremely contaminated, Zn was strong to extremely contaminated and Cu was strongly contaminated in their geoaccumulation study. Moreover, Radomirović et al. 75  www.nature.com/scientificreports/ polluted soil by Zn based on their geoaccumulation study at a former painting industry facility. Furthermore, the PI and PLI depict that all the samples (both levels) are unpolluted but surface soil is more polluted than the sub-surface. Hu et al. 11 also found a safe level of HM pollution in the soil at a coastal industrial city in the Yangtze River Delta, China. The soils have a low level of contamination due to the lower concentration of HM in collected soils. The distribution of HM and their pollution level depend on natural and anthropogenic factors such as rock weathering, urban-industrial wastewater, transportation, and modern agriculture 4,13,33 . Moreover, Wang and Zhang 12 found that the deposition of HM in urban roadside soil was due to restoring damaged roads and maintaining green belts. However, the present study found that the distribution of HM and their pollution level are influenced by modern agricultural practices, soil structure, precipitation, and river water fluctuations. The HM pollution concentration is lower in the soil profile because of two factors-(1) sandy soil of the char and (2) frequent char inundation. The sandy soils absorb the pollutants and infiltrating water through the soil layers leaches pollutants down to the deeper soil layer or the river (Fig. 11d,e). A similar finding was reported by   76 . Moreover, the low-lying parts of the chars frequently become inundated by monsoon floods almost every year that disturb the geo-accumulation process resulting in a lesser pollution concentration. Our study findings are also supported by Sarkar et al. 52 . HM pollution of the soil though low is contributed mainly by the charland agricultural practice pre-monsoon season because of the huge application of chemical fertilizers and pesticides for agricultural production (Fig. 11f).

found moderately
The study reveals that the soils have a low ecological risk and low potentiality of ecological risk. However, the ecological risk is comparatively higher in the surface soil than in the sub-surface soil due to the lower concentration of HM in the sandy soils. Additionally, the potential ecological risk of HM is ranked as Cu > Zn > Mn for both levels of the soil. Liu et al. 77 found that most of the HM (As, Cr, Cu, Ni and Zn) had a low potential ecological risk in coal mining areas in China. Additionally, Perumal et al. 33 also found a low-potential ecological risk from their studied HM and considerable ecological risk from Cu and Zn, and similar results were also found by Sarma Bora et al. 74 . Besides, Radomirović et al. 75 found moderate pollution and ecological risk for most soils at a former painting industry facility. Moreover, as per the enrichment factor, the HM is ranked as Mn > Zn > Cu for both levels of soil. The HM enrichment in sub-surface soil ranges from very severe to extremely severe while they range from severe to extremely severe for the surface soils. Perumal et al. 33 found no enrichment for Mn, Fe, and Cu but Radomirović et al. 75 found a moderate soil enrichment for Zn in the soil. Sarma Bora et al. 74 found significant enrichment for Cu and Zn in solid waste dumping site soil near the Morabharal River, Tezpur town, India.
Moreover, the final relative closeness score of the ideal solution indicates that the sub-surface soils have lower soil pollution than the surface soils in the present investigation. Further from the geostatistical modeling, it was observed that the SK technique was expected to be the most accurate interpolation. Islam et al. 18 found that the SK technique was the most accurate interpolation model for As, Mn, Zn, and AI while the OK technique for Fe and IDW with power1 for Ba. The kriging models (SK and OK) found as the most accurate interpolation technique for the spatial distribution map of all elements. Islam et al. 18 also found a similar result for geostatistical modeling analysis in their study in the Rangpur district, Bangladesh. Moreover, the Gaussian semivariogram model was observed as the best-fit semivariogram model among the other semivariogram models (circular, spherical and exponential) studied.
Hence, this investigation would be helpful for decision-making for sustainable agricultural development, improving the degrading riverine ecology and better representation and map-making processes development. Although this present study has some limitations. The main limitation of the study is data unavailability of other trace elements or HM like arsenic (As), nickel (Ni), cadmium (Cd), molybdenum (Mo), chromium (Cr), mercury (Hg), lead (Pb), fluorine (F), and selenium (Se) and lack of sufficient fund for testing these trace elements. Despite those drawbacks, the finding of this study may be helpful to stakeholders including regional planners, and agricultural engineers in the DRB.

Conclusions
The results of this study portray the nature of PC parameters, HM pollution, and ecological risk in mid-channel bar surface and sub-surface soils from the DRB. This study also shows the best-fit spatial distribution models for the studied variables from both levels of soil. The major findings are mentioned below.
• The concentration of mean pH shows that surface soils are moderately acidic and the sub-surface soils are slightly acidic. Moreover, the Ficklin diagram shows that 76.67% of surface soil is acid-high metal while 83.33% of sub-surface soils are near natural high metal. • Shepard Triangle Diagram shows that 83.34% of surface soils and 90% of sub-surface soils are sandy loam in the mid-channel bar of Damodar River. • Similar results were observed as the HM concentration in soils was low, the contamination level was also low (based on the CF, Cd, and Cp values) and the PI and PLI also indicate that all the samples (both levels) were the unpolluted or low level of pollution but the geoaccumulation index shows that all the soils are uncontaminated to moderately contaminate; • The potential ecological risk index also shows that all the soils (both levels) have a low level of ecological risk.
Inverse results were observed between potential ecological risk and enrichment as the potential ecological risk of HM is ranked as Cu > Zn > Mn while enrichment as Mn > Zn > Cu for both levels of soils. • Moreover, TOPSIS shows the sub-surface soils have lower soil pollution than the surface soils of the studied chars. • Geostatistical modeling analysis shows that the SK technique was expected as the most accurate interpolation model as 61.54% and 76.92% of the surface and sub-surface soil sample's variables fit with the SK interpolation technique compared to other interpolation techniques (OK and IDW) and Gaussian semi-variogram model was the best-fit semi-variogram model among the other semivariogram models (circular, spherical and exponential) studied. • The results of numerous indices show that the pollution or contamination level and ecological risk from surface soil are more as compared to the sub-surface soil while the enrichment, geoacumulation, and potential contamination of sub-surface soils HM are high than the surface. Moreover, the sub-surface soil is more sandy than the surface soil and the concentration of Cu and pH is high than the surface soil. The SK technique and Gaussian semi-variogram model was expected to better fit the model as the most accurate interpolation model for sub-surface soil variables compared to the surface soil variables.